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A bs  tract 

After  developing  the  formal  integration  of  S IB /at  = s?/  xB) , 
it  is  3hown  that  cyclonic  convective  motions  in  the  core  oroduce 
magnetic  loons  in  meridional  planes  through  interaction  with  the 
toroidal  magnetic  fields  Expressing  these  loons  in  terms  of  the 
usual  orthogonal  vector  modes,  it  is  shown  that  they  result  in  a 
Dredominantly  dipole  field.  Together  with  the  nonuniform  rotation 
of  the  core,  which  produces  the  toroidal  field  from  the  dipole 
field,  the  cyclonic  motions  result  in  a comnlete  s elf-regenera- 
ting magnetic  dvnamo.  We  conclude  that  anv  rotating,  convectinp, 
electrically  conducting  bodv  of  sufficient  size  will  nossess  a 
magnetic  field  generated  by  this  dvnamo  mechanism.  The  nossibility 
of  an  a brunt  reversal  of  the  field  is  discussed. 
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1.  Introduction 

In  ths  magnetohydrodynamic  theory  of  the  Earth’s 
magnetic  field  (Elsasser,  1946j  Bullard,  1949),  it  is  shown 
that  the  nonuniform  rotation  of  the  liquid  core  of  the  Earth  will 
'produce  a strong  toroidal  field,  (that  is,  a field  along  the 
circles  of  latitude).  The  nonuniform  rotation  of  the  core  is  the 
result  of  the  Coriolis  forces  which  act  upon  the  rising  and  falling 
convective  motions.  If  we  are  to  show  that  the  dipole  field  is 
due  to  a self-sustaining  dynamo,  we  must  demonstrate  the  existence 
of  a feedback  link,  i.e.  of  the  mechanism  by  which  the  dipole  field 
is  regenerated  from  the  toroidal  field. 

Cowling  (1953)  has  shown  that  there  is  no  very  direct 
way  of  producing  e dipole  field  with  simple  fluid  motion}  we 
must  in  particular  abandon  all  models  possessing  rotational 
symmetry.  Cowling’ 3 point  is  that  in  a dipole  field  there  must 
be  a line  singularity  about  which  the  magnetic  field  circles. 

The  field  vanishes  at  the  singularity,  but  the  curl  does  not. 

But  for  a steady  state  field,  the  field  equations  tell  us  that 
the  field  IB,  curl  B , and  the  current  density  ji  , must  all  become 
small  of  the  same  order  if  one  of  these  quantities  is  small. 

Hence  such  a field  cannot  be  maintained  in  a steady  state.  The 
linearity  in  B implies,  furthermore,  that  a time  dependent 
velocity  field  is  equally  unable  to  maintain  or  amplify  the  field 
in  the  average. 

In  this  paper  we  first  point  out  that,  owing  to  the 
Coriolis  force,  the  radial  convective  streams  must  produce  local 
cyclonic  and  anticyolonic  circulations  somewhat  resembling  those 


observed  in  the  atmosphere.  The  deformation  of  the  toroidal  field 
by  these  motions  can  be  shown  to  give  rise  to  the  formation  of 
loops  of  the  magnetic  lines  of  force  in  the  meridional  planes  which 
coalesce  to  a moan  meridional  field.  The  dipole  can  be  shown  to 
be  by  far  the  largest  harmonic  component  of  this  field.  Thus  the 
interaction  of  the  cyclonic  and  anticyclonic  local  motions  with 
the  toroidal  field  provides  the  desired  feedback  linkj  together 
with  the  interaction  of  the  nonuniform  rotation  of  the  core  and 
the  dipole  field,  this  constitutes  a complete  regenerative  cycle. 
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2.  Dynamics 

The  primary  motions  in  the  core  are  assumed  to  be 
convective  motions  in  a radial  direction,  (Bullard,  1949$ 

Elsasser,  1950)  r It  is  readily  demonstrated  that  they  result  in 
a nonuniform  rotation  of  the  core- in  which  the  outer  layers  of  the 
core  lag  behind  the  average  rotation  and  the  inner  ones  exceed 
the  average. 

In  the  atmosphere  of  the  Eahth  the  Coriolis  force  of 
the  large  scale  cyclonic  motions  shears  the  vertical  motion  to 
such  an  extent  that  the  actual  direction  of  flew  in  a rising  eddy 
may  be  less  than  a degree  from  the  horizontal.  We  do  not  expect 
the  shearing  to  be  so  large  in  the  core  because  of  the  retarding 
forces  exerted  by  the  toroidal  field.  Bullard  (1954)  has  shown 
that  for  any  self-sustaining  dynamo  there  exists  a steady  state 
velocity  which,  when  exceeded,  causes  the  magnetic  field  to  grow 
until  its  reactive  forces  balance  the  mechanical  forces  driving 
the  dynamo.  Assuming  that  the  core  of  the  Earth  Is  a self- 
sustaining  dynamo,  we  expect  to  find  It  operating  near  the  steady 
state  velocity  with  the  Coriolis  force  balancing  the  reaction 
of  the  toroidal  field. 

Defining  the  magnetic  Reynolds  number  if  the  nonuniform 
rotation  In  terms  of  the  radius  of  the  core,  Bullard  (1954)  finds 
it  to  be  of  the  order  of  25  for  steady  state  operation.  If  we 
choose  6 = 3.5  x 104  mho/m  (Elsasser,  1950)  we  obtain  toroidal 
velocities  of  the  order  of  0.1  ram/sec.  T0  obtain  an  estirrate  of 
the  radial  velocities  required  we  shall  assume  a suitable  value 
for  the  toroidal  field  and  estimate  the  ^oriolis  force  needed  to 


overcome  its  reaction.  A value  of  40  gauss  seems  a reasonable 
lower  limit  on  the  toroidal  field,  , and  four  gauss  for  the 
dipole  field  Bj  . The  reactive  force  per  unit  mas3  is  (V  * B)x 
newtons/kg  where  is  the  density  of  the  medium.  Approximating 

this  as  By  where  L is  some  characteristic  length,  we 

obtaiin  a force  of  the  order  of  1.3  x 10-^  knewtons/kg  for 
/o  - 1C  gm/cm^  and  L = 1000  km.  With  ' — - 10-4  sec  ± we  obtain 

a lower  limit  on  the  average  radial  velocity  of  10  mm/sec. 

The  individual  active  convective  regions  probably  have  radial 

_2 

velocities  one  or  two  powers  of  ten  larger  than  this,  say  10 
to  10_1  mm/sec,  and  we  estimate  the  ratio  of  toroidal  to  radial 
velocity  to  be  ten  or  less,  whereas  in  the  atmosphere  this  ratio 
is  generally  of  the  order  of  a hundred. 

In  this  paper  we  shall  disregard  the  nonuniform  rotation 
and  shall  consider  only  the  interaction  of  the  radial  motions  with 
the  toroidal  field.  Consider,  then,  a radial,  convective  stream 
wioh  the  associated  influx  and  efflux  of  matter  at  its  ends,  A 
rising  stream  will  result  in  an  efflux  of  fluid  near  the  surface 
of  the  core  and  an  influx  near  the  center,  and  conversely  for  a 
descending  stream.  As  is  well  known  from  meteorological  phenomena, 
the  Coriolis  forces  resulting  from  the  influx  and  efflux  produce 
a rotation  of  the  convective  column  about  its  axis.  T0  formulate 
this,  we  define  a local  cartesian  coordinate  system  (§ , ) 

where  the  f -axis  is  tangent  to  the  line  of  longitude  through 
the  origin  and  the  positive  £ -axis  points  south;  the  ^ -axis 
is  tangent  to  a line  of  latitude  through  the  origin  and  points 
east,  and  the  jf  -axis  ’ ~ directed  radially  outwards.  We  define 
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and  the  pcuar  angle  ^ measuring  azimuth  about  the  jf  -axis  from 
the  § -axis.  In  this  system  we  represent  an  outflow  from  the 
radial  motion  as 

v,  = v,  (_  cos  +-  -sir  yW  } 5 1» ^ ^ — ,v//<>  j '°5^  1 S 

where  the  € ; s are  unit  vectors  in  the  corresponding  directions. 

The  angular  velocity  of  the  Earth  is 

<•**  — c^_)  ( CosS  --  C.  sm  ©^ 

where  9 is  the  colatitude  of  the  origin  of  the  ( i, n , / } system. 

The  Coriolis  force  due  to  this  efflux  is 

— 2 C XV,j  = 2 U\',(.*)|  «>  CoS©  CO  ?></'  ~ C,  CO'S©  Cos  yl  - Cj,  -S/o  © S>iH 

The  y ^component  of  the  curl  of  this  force  is 

cur  |y  \-Z(  a***,)).  = JL  (- 2 & x yg,\ 

= * 2 w CoS  © — £/»  U,  ) 

)/J^  Is  positive  for  an  efflux  and  negative  for  an  influx} 
this  defines  the  sense  of  the  circulation. 

We  shall  furthermore  assume  that  the  dissipative  forces 
of  viscosity,,  turbulence,  etc.  may  be  neglected  in  a first  approxi- 
mation (the  Reynolds  number  being  large)  with  the  result  that  the 
convecting  regions  will  always  show  a sense  of  circulation  corres- 
ponding to  an  influx i The  convective  region  is  set  in  rotation  by 
the  Influx}  In  the  absence  of  dissipation  the  angular  momentum  of 
tne  convective  column  Is  conserved  and  is  stopped  only  by  the 
equal  and  opposite  angular  Impulse  of  the  efflux  at  the  terminus 
of  the  convective  stream.  In  the  northern  hemisphere  the  circula- 
tion Is  counterclockwise,,  in  the  southern  hemisphere  clockwise.,  if 
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viewed  from  outside  the  core.  These  motions  are  illustrated  in 
figure  i. 


3.  Magnetic  Field 


Bullard  (1954)  has  analyzed  the  toroidal  magnetic 
field  generated  by  the  interaction  of  the  nonuniform  rotation  with 
the  dipole  field.  The  result  essential  to  the  feedback  link 
discussed  in  this  paper  is  that  the  toroidal  field  has  a maximum 
about  half  or  two  thirds  of  the  way  out  from  the  center  of  the 
core  and  drops  off  to  zero  near  the  center  and  near-  the  surface. 

The  conductivity  of  the  core  is  sufficiently  high  so  that 
fluid  velocities  of  10  “ mm/sec  or  more  may  be  expected  to  carry 
the  lines  of  force  along  bodily.  Thus,  a conve  ctive  upwelling 
produces  an  upward  bulge  in  the  toroidal  field.  The  rotation  of 
the  upwelling  about  the  radial  direction  twists  this  bulge  into 
a loop  with  a nonvanishing  projection  on  the  meridional  plane 
resulting,  in  the  northern  hemisphere,  in  the  sequence  illustrated 
in  figure  2.  Figure  3 shows  the  similar  process  occurring  in  a 
sinking  column  of  fluid  in  the  northern  hemisphere.  In  the  southern 
hemisphere  the  se:  ie  of  the  toroidal  field  arid  of  the  rotation  of 
the  convective  column  are  reversed.  4 illustrates  the 

complete  dynamo  model;  A typical  line  of  force  of  the  dipole  field 
of  the  Earth  is  shown  as  having  been  drawn  out  by  the  nonuniform 
rotation  to  give  a toroidal  field}  the  loops  produced  from  the 
toroidal  field  by  convective  motions  in  both  hemispheres  are 
shown.  We  note  that  the  loops  produced  nearest  the  surface  of 
the  core  arise  from  rising  convective  currents  and  have  the  same 
sense  of  circulation  as  the  overall  external  dipole  field  of  the 
core}  the  loops  produced  near  the  center  of  the  core  have  the 


opposite  sense,  but  their  great  depth  means  that  they  will  decay 
before  diffusing  to  the  upper  part  of  the  core  where  the  main 
amplification  of  the  toroidal  field  occurs.  The  dipole  field 
observed  outside  the  core,  which  is  due  to  the  coalescehca  of 
all  the  loops  formed  by  the  convective  motions,  will  have  the  same 
sense  as  the  original  dipole  field. 

This  qualitative  discussion  indicates  that  the  simple 
dynamical  model  of  the  convective  motions  presented  in  the  previous 
section  should  result  in  a self-sustaining  dynamo  with  a dipole 
field  . Insofar  as  our  model  is  correct, it  would  follow  that  any 
suf f ic iently  large,  convecting,  electrically  conducting  fluid  body 
should  exhibit  a magnetic  dipole  roughly  parallel  to  the  axis  of 
rotation  if  the  body  rotates  rapidly  enough.  On  reversing  the 
sense  of  either  the  magnetic  field  or  the  rotation  of  the  body 
in  the  diagrams  we  see  that  the  amplif icatory  process  is  inde- 
pendent of  either.  Thus,  the  sense  of  the  dipole  field  will 
depend  only  on  the  sense  of  the  initial  field  from  which  it  was 
generated* 


The  entire  process  has  been  dubbed  the  bathtub  effect 
for,  as  every  physics  freshman  has  heard,  the  rotation  of  the 
water  running  into  the  drain  of  a bathtub  Is  the  result  of  the 
Coriolis  force. 


4.  Generation  of  Loops 


The  basic  equations  for  the  deformation  of  the  magnetic 
field  JB  by  the  velocity  field  \g  are,  for  a perfectly  conducting 
fluid  (Elsasser,  1950) 

= Vx(^xB)  U) 

o t ~~ 

Thi3  diff ex-ential  equation  is  equivalent  to  the  integral  equation 


at  [b-js^o 

where  dS  is  an  element  of  a surface  moving  with  the  fluid, 
fellows  that  the  complete  integral  of  (1)  is 


/ 


B-dS 


cons 


f&ni 


(2) 

It 

(3) 


which  may  be  written  as  three  linear  simultaneous  algebraic 
equations  for  B;  by  choosing  infinitesimal  areas  over  which  to 
Integrate.  Thus  we  consider  the  infinitesimal  element  of  volume 
determined  at  time  t = 0 by  the  three  linearly  independents  but 
otherwise  arbitrary,  vectors  6a;  J <5>bt-  £ct  . Let  now  one  corner 
of  this  element  of  volume  be  at  , then 


ba ; = x,  * — x 
At  some  later  time. 


‘ . — x.'  , & = y.,  ‘ - ; 

t , we  write  for  the  same  material  element 


&a.  - x. ' - X.'  , &B,  -X.‘ -x:  , %c,  = A,'  - x.‘ 

The  Lagrangian  coordinates  X are  functions  of  the  initial 
position  and  of  t,  so  that 

X*~  X‘C *\0 


(5) 


Now 


I V 

i 
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X.'  - XV.O  -XV*  - xWfr-  ^ 


s&»  +-  Ox(sai'- 


and  similarly  for  X,  and  Xa  • Hence,  by  (5) 

feA,  = SaK;%Bt  = y SO  - ^4  ( 6 ) 

The  areas  of  the  faces  of  the  element  of  volume  are 
(we  use  summation  convention) 

where  t‘,k  is  the  usual  per*mutation  tensor,  being  +1  if  i. j k is 
an  even,  -1  if  it  is  an  odd  permutation  of  1,  2,  3,  and  zero  other- 
wise. If  the  magnetic  field  is  written  as  (X\  t)  with  the 
initial  value  B>;  o)  » b ■ , the  flux  through  the 

face  bounded  by  &A;  and  “SB;  is  initially  b>c  & A . SBx  5 at 

time  t it  Is  B;  £.  & A\  SB.  1 and  similarly  for  the  other 

twn  aides.  Thus,  upon  using  (6)  to  eliminate  SBi  and  j 

we  find  from  (3) 


(7) 


b;  <£.  ‘iU  &*,  S b*  = 

B;  «£:  sa, 

■£>*.-  2>x.* 

%»b;  ^C|  - 

g . %\or 

be*'1*  €>c,  % a * ~ 

ts . e sL^L  g-X  ^Cr 

' -fcx* 

Since  the  &j;  ; Sio;  and 

Sec  are  essentially 

may  set 

*•  S,;  i S b ,• 

- ; S Ci  = Sti 

where  SCj  Is  the  usual  Kronecker  delta.  Then  (7)  reduces  to 

u _ R &2L'  b = B;  , fc>.  = Bl  (8) 

These  three  linear  algebraic  equations  can  be  solved 
to  yield  B,  (X  ,i)  in  terms  of  and  the  Lagrangian 
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I 

I 

i 


i 

I 

1 


coordinates  X . To  facilitate  this  solution  we  denote  by  J,;  the 
derivative  X $ the  Jtj  form  a matrix  which  we  shall  write 

as  (J  ).  The  determinant  of  this  matrix,  denoted  by  IJ!  , is  the 
Jacobian  of  the  transformation  X = X ‘ (*■  O * If  ( b ) and  ( B) 
denote  the  column  matrices  with  components  b,  and  B,  respectively, 
then  (28)  may  be  written 

(b)  «=  J)  (B) 


where  the  adjoint  matrix  is  defined  as  usual,  by  the  relation 

*jjj  = (j;"  iji 


and  hence 


Thus,  finally 


1) 


(b;  = 


B;  = 


) T I 


U )ct 


b ■ 


O) 


(10) 


\J\  ax- 

which  is  the  general  solution  of  (1). 

We  shall  now  use  (10)  to  demonstrate  that  rotating 
convective  currents  will  produce  loops  of  magnetic  flux  which  are 
displaced,  relative  to  the  height  of  the  toroidal  field,  in  the 
direction  of  the  convective  motion.  Since  we  use  the  I-agrangian 
method,  the  velocity  field  must  be  expressed  in  terms  of  the  dis- 
placement of  the  individual  fluid  particle  s.  On  the  assumptions 
made  the  trajectory  of  a particle  will,  in  the  first  approximation, 
be  a helix.  We  set 


dMt r,  C,, 

dt  a x - 


v ' 1 ~7t  ~ "a  s S(^J  , gf  = 

The  pitch  of  each  helix  is  (y./ w-X Rq^)/  S(^>) 


(11) 


Tj  Some  steps  of  this  integration  procedure  have  already 
been  indicated  by  Lundquist  (1952) 


-J 
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Integrating  (11)  we  obtain 

= X , i - A c„[ff  5W1  * xj  A«  A )ff  50.  >U  X.]  . / - w ♦ X , 

where  \ , > and  Aj  are  constants  of  integration.  If  § 

and  J'  = J£  for  t = 0,  then 


X 

-VC  €.: 

■*VJ  , 

X.  = _ 

L , c 

ot 

A = 

' ■ a 

€. 

1- 

It  f 

ollows 

that 

II 

II 

> • 

X tr' 
m ft 

VX 

x > 2>  r. 

= o J 

~sX. 

SB 

h X, 

= o 

> 

“ci  X !. 

_ *v . 

& x. 

— j « 

■sX. 

bf. 

XJ  1 

" V 

> 

2>J. 

— O 

A, 

■2>jr. 


Using  (14)  it  is  readily  shown  from  (12)  that 

i§l  = +■  wJ  - x §o  s'(a)  t 4 


If.*  i C\.5  ^ v s*(»t^ 


A 

. '-b-|  = o 

' &jr. 


(12) 

(13) 


(14) 


li.  = " i - £ O + f-  «. 


1^.=  ^ M *,Vr0  + ^ 'x.'Z'WlJL  , 11a  = 

-3. 


o 


- v.  f?'OM  i , 2.J 

**-  > ' -2.1..  A 


The  flow  is  incompressible  and  ! J - i . 

We  now  specify  our  magnetic  field  as  purely  toroidal 
and  locally  a function  of  height  only*  say 

(15) 


Bf  = By  = ° , B,=  B(jr; 


Prom  (10)  we  then  find 


"13" 


b,  - B«)Ji  = b(x)[-i  (vi  - f.  x)-  f v sou  f ] 

b,  = = B(rjf-7:(ts  t-  f*  s-(m  i7  <16> 

"b^..  *-  ^ A J 

By  = - B(rjv.  R'wt  4^ 


Let  us  now  use  (12)  to  eliminate  ^ ^ and  X from 

(16).  From  (13)  wa  may  rewrite  (12)  as 


§ — § c < 


> S(^‘]  - 1.  v.f^  s(\it] 


a = s. s 


in[-£scx)t ] + i\.  co';‘  [ I‘  5(A)  t] 


Solving  these  two  equations  for  §.  and  and  using  (12)  again 


we  obtain 


~ £ cos  [£  bcah]  sca;(] 


1.  — _ S sil 


[2-  Swi]  SOOl] 


X - Jr  - * roo  t 

We  let  t = Xcos>A  and  r^=  A am  ^ so  that 


i.  = Uos^-^  soot]  ) >\.  = a SGOt] 

Using  (IV)  and  (18)  we  now  rewrite  (16)  as 

= -Sa-v..F?(A)t){sin(^.  SWJ)  + f*  S(>J  >t  srn^-  ^SCA)4]} 


bv  =+b(x-v.  R,(A)t)|cos^s(A)t)  + ^ s'(a)au6j^  a-.  sa;tj^19) 

5y  =+B  (Jf  - X'.ROJt)  v.  R,'(>)t  sm[^-  f^SC X)t]  J 


-*I4 


Let  us  investigate  the  field  in  tie  jff  --  plane.  We 


put  = o and  obtain 
Bf  = -B(jT-v.  r?(»0  -&sc»t] 


-v.  R(» 


I 

t){c°s[^s(A)t]  - f*  j>  (20) 


«-  e(*  -v.R(x)t)  x t si»-,  f*.  scx)t]  y 

»Ve  are  primarily  interested  i:n  the  loops  formed  and  we 
therefore  subtract  from  this  field  the  portions  that  do  not  con- 
tribute to  the  loops,,  namely  r those  due  to  the  motion  in  the  Jf 
direction  alone  and  those  due  to  the  circular  motion  about  the 
y — axis  alone,,  so  as  to  retain  only  the  part  that  results  from 
the  superposition  of  these  motions.  Thus,,  for  w.  = o we  write  the 


'ic  Id  as 


<*1  = O , ^ ; a;=o 


For  v.  = o we  write 


/Gf  = - B(  X)  si*  sc.ot] 


A - ^ B(r3  [co5[^  5(A;tJ  - ^ S(AJ  M sirvf^i  S(X)1 


A = O 


Hence 


Bf-  = [B(jf3-  B(j  _ v.  siefs*  sa)t] 

B,-^-A=  fB0-v.R(X)O-  BUj{,os[BS0>)t}-£  %-MU  Si„[?soo4| 

- B(J  - v.  fR(X)0 

B,  " **  "A  = - B (jr  - v.  R(x:t)v.  RVx.T  s/nt] 

Z.  - — ' • J 


I 


The  projections  of  the  lines  of  force  on  the  $ § -plane 
satisfy  the  differential  equation 


djr  _ Br  - = 


Bff  - *:  -ft 


BCJr-v.  R(*)i)  v.  R'mt 
2>  ( jf  - v.  R(X)t)  - S>  (J) 


We  define  the  variable  ^u.  as 

= v.  RCA)  t 


so  that 


Then 


«■  v.  fR(A)t*JA 


J5 


This  nay  be  rewritten  as 


We  define  the  function  X(x)  as 

X (.x)  — B 

The  integral  of  (23)  becomes 

1(5-  v.  *C*)0  *I»)tC 


(21) 


(22) 


( 23) 


(24) 


(25) 


If  we  assume  that  RCA)  is  a function  having  a maximum 
at  A = 0 and  decreasing  monotonically  to  zero  for  large  | A|  ,11 
can  be  shown  that  (25)  results  in  magnetic  loons ; the  loo os  for 
special  case  that 

F?o>)=  e*p(-  y/a*)  , e>cs)=  s/*') 

are  shown  in  figure  5 for  v.  t = a . The  neutral  noint  is 


%=  O J X = v.t/2  = a/2. 


-*• 
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Thus,  we  have  demonstrated  thut  the  loons  are  cisnl seed 
in  the  direction  of  the  convective  motion. 

It  is  well  to  note  that  if  v-/.  is  so  large  that  the 
convective  column  rotates  more  than  one  half  a revolution  after 
passing  through  the  toroidal  field,  the  sense  of  part  of  each  loon 
will  be  reversed.  We  shall  assume  that  loops  in  the  first  half 
revolution  predominate  and  the  average  sense  is  as  illustrated  in 
figures  2.  3,  and  4. 

In  the  next  section  where  we  study  the  coalescence  of 
loons;  we  shall  see  that  a set  of  loops  of  one  sense  of  flux  above 
the  toroidal  field  and  another  set  of  the  opposite  sense  below  this 
field  combine  in  such  a wav  that  the  net  average  field  is  dominated 
by  the  loops  at  the  higher  level,  i.e.,  the  resulting  noloidal  field 
of  the  core  has  the  same  sense  as  the  upper  "'cops.  We  know  little 
about  the  dynamics  of  convective  motion  and  it  is  quite  possible 
that  in  the  core  there  are  rising  and  descending  streams  in  equal 
numbers  and  of  symmetrical  structure.  Again,  it  is  possible  that 
the  convecti  VC  H30 t ion  is  assymmetrical , the  ascending  streams 
being  concentrated  in  narrow  regions  and  the  compensating  descending 
motions  being  spread  over  a large  area  (the  familiar  pattern  of 
atmospheric  thunder  storms).  Now  it  can  he  demonstrated  that  in 
such  a model  of  convection  only  the  intense  rising  currents  produce 
oi-^redabl  e magnetic  Iccr.s  whereas  the  contribution  of  the  spread- 
out  subsiding  motions  to  the  meridional  ma /netic  field  is  negligible. 
This  might  be  shown  on  using  the  preceding  analysis,  but  the  cal- 
culations turn  out  to  be  extremely  involved.  In  the  appendix 
we  therefore  solve  equation  (l)  again,  this  time  by  a perturbation 
procedure.  In  the  approximation  which  we  use;  the  displacement  of 
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the  loops  in  the  direction  of  the  convective  str-f  am  does  not  vet 
anear^  on  the  other  hand  we  obtain  simple  analytical  expressions 
for  the  secondary  magnetic  field,  and  these  exhibit  clearly  the 
fact  that  the  contribution  of  the  snread-out  descending  motions  to 
the  meridional  loons  is  negligible-  The  general  conclusion  is  that 
we  can  construct  a self-sustaining  dynamo  on  the  basis  of  either  a 
svmmetrical  or  an  asymmetrical  kinematical  model  of  convection. 
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5.  Coalescence  of  Loops 

It  remains  to  be  shown  that  a number*  of  loops  with  the 
sense  indicated  in  figure  4 will  actually  result  in  a field  ob- 
served outside  the  core  which  is  predominantly  a dipole  field. 

Tne  study  of  the  geomagnetic  secular  variation  -shows  (Elsasser, 

1950)  that  there  are  about  fourteen  distinct  local  irregular! t ies 
in  the  field  of  the  Earth.  Taking  this  figure  as  a lower  limit  on 
the  total  number  of  loons  in  various  stages  of  decay  in  the  core, 
the  total  number  of  Ioodb  might  be  as  high  as  twice  t-.his  amount. 

We  shall  assume  that  the  loops  are  distributed  randomly  over  a 
sphere  of  given  radius  concentric  with  the  core,  and  we  shall 
represent  this  distribution  by  a uniform  average  density. 

Before  we  carry  out  the  calculations,  let  us  briefly  dis- 
cuss this  averaging  nrocess.  We  shall  find  that  the  looos  produced 
beneath  the  toroidal  f ield  by  descending  currents  of  fluid  (and 
which  have  the  ’’wrong”  sense  for  feedback)  are  buried  so  deeply 
in  the  core  that  they  do  not  contribute  appreciably  to  the  resultant 
mean  poloidal  field.  Since  the  coalescence  problem  is  linear,  we 
may  demonstrate  this  fact  in  the  following  fashion:  We  first  com- 

pute the  magnitude  of  the  mean  noloidal  field  with  the  restriction 
that  the  looos  are  confined  within  a sphere  of  arbitrary  radius 
not  exceeding  the  radius  of  the  core.  We  can  then  represent  the 
condition  where  there  are  loops  of  one  sense  in  the  upper  and  loops 
of  the  opposite  sense  in  the  lower  part  of  the  core  by  means  of 
a linear  superposition  of  two  concentric  spheres  of  different 
diameters,  each  being  filled  with  loops  of  only  one  sign. 
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Let  cf*  be  a system  of  polar  coordinates  whose  oi-igin  is  at 

the  center  of  the  core*  Let  a designate  the  linear  dimensions  of 
a tynical  loopj  we  assume  the  loons  small,  a«  P where  P designates 
now  the  radius  of  the  core.  Consider  a loop  whose  center  is  at 
r=  Ft*,2>*L=  cS , * = <^3  For  any  point  on  the  loon  we  set 

r = -iyO  } 2>  * = gca-i  q>  * =.  -h  cf  ( 26  ) 

where  /a,  ib,  are  small  quantities  of  order  a/fy1  whose  higher 

f 

nowers  will  be  neglected. 

To  simnlify  the  calculations  we  next  assume  that  the  loop 
can  be  described  by  means  of  a gaussian  svmnetrical  about  the 
loop  center.  We  chooso 

B_  = — Ei.  Jp  7>  exp  f~i  (/*  '♦  R%i>'  + R' tfi1 jj 

= +■  B.  ^ f si^v  e>  v 

B,  = o . 

the  lines  of  force  are  circles  with  centers  on  the  cp-axis;  together 
they  form  a torus  of  magnetic  flux.  Some  of  the  later  analysis 
is  more  conveniently  carried  out  if  we  work  with  the  curl  of  the 
magnetic  field,  since  this  is  a toroidal  rather  than  a noloidal 
vector  (Elcasser  1946).  We  shall  need  onlv  the  (p  -component  of 
the  curl  for  reasons  that  will  annear  presently. 

(VxbJ^  = (/—  ^‘+  RV'  t-  R>*si«'05  cW] 

V/e  now  expand  (27)  in  terms  of  orthogonal  toroidal  vector 
modes,  v.'horeupon  we  average  over  all  possible  positions  of  the  loop 
centers.  If  the  distribution  function  of  the  loons  is  independent 
of  <£>  g only  the  rotationally  symmetric  modes  (zonal  harmonics) 
will  survive  the  averaging  process.  The  toroidal  zonal  modes  are 
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* *-V 


where  j„  (.x)  is  a spherical  3ossel  function  defined  (Stratton,  1941) 
as  (n/^Z  *-))  ("*■)  and  where  the  boundary  conditions  require  the 

to  be  the  roots  of  (Elsasser,  1946) 

(k..F»)«oj  <29) 

P beinp,  the  radius  of  the  core. 

We  now  develop  (28)  in  ascending  rowers  of  the  small  quan- 
tities and  ~i}  defined  in  (26). 

"T^u  = (Am+  ES,.^  + C.  > + +o(^A)  (30) 

We  will  be  interested  only  in  D„,  and  FT,  t which  are  readily  found 
from  (28)  to  be 

Ft . « -L  Jp  (<os<ep)  FT.  « g;  (31) 


JR' 

The  modes 


TS * 

are  orthogonal  so  that  if  we  write 


V x B = T 


multiply  by  X...  and  integrate  over  the  volume  of  the  core,  we 


obtain 


fff**  ^S)-X,  = Xf  X., 

Noting  that 

T.  * ' 

i j v.**  r j __  f , /,  ^ - 1 1 2^ 

Jdv  SI^  [_  — pr- j - R c-)J  = —^7- 

(MacRobert,  1948,  pg.  105),  and  that 

-*  = t /dr  = i p'j.(k.pj 

(from  (29)  and  Jahnke-iinde , 1945,  pg.  146),  we  obtain 


X-X.,  “ PJi,'(  tvP)  * 


From  (27)  and  (3C)  we  find 
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j y^dJr  (5?/  E>)  • 


T„>  = - B 


. cM  a « p 


D + ^ t O 
M + R 


(l%)i 


(34) 


all  other  terms  cancelling  out  for  reasons  of  symmetry.  Sub- 
stituting (33)  and  (34)  into  (32)  we  obtain  finally 

FI.  1 a4  1 ) (35) 


- - B. 


(*• 


R**/  2n(n+i) 

We  next  consider  the  tine  dependence  of  the  loops.  They 

will  ancear  in  some  random  fashion  and  decay  subsequently.  The 

easiest  way  to  discuss  the  decay  of  the  normal  modes  will,  be  to 

use  the  time  dependent  differential  equation  with  a suitable 

source  function  fn,  C^)  • The  amplitude  of  the  mode 

then  satisfies 


A + -txK.U)  = L (t)  , A,.  = jXTcT 


(36) 


The  general  solution  of  this  equation  is 

?r.(i/  = JcJt'  exp  0] 

Let  us  assume  that  the  loops  anoear  suddenly  at  random  times 
tf  3 thereafter  they  merely  spread  out  by  diffusion.  It  is  readily 
seen  that  the  source  function  T„t({)  is  nrooorti  onal  to  the  rate 
of  appearance  of  the  loops j we  therefore  take 


k*  = c..  y~  & ^t-  to 

where  S>Ct~t,0  is  a -Dirac  delta  function.  Then 

?r.ct)  = y~  [a,.  (t:  ~t)J 

where  the  summation  is  over  all  i for  which  t > t; 
The  average  of  ^ ( t)  over  time  is 

* T 

<*■>  - rl '<«) 

which  gives  after  some  calculation 


/ 


im 


I 


. S _!_  — 

— A^s  A», 
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(37) 


T— » co  ^ T 

where  u is  the  number  of  loons  nroduced  ner  un^'t  time. 

We  assume  the  apriori  probability  of  a loop  appearing  at 
to  be  K (F?)  , normalized  to  rive  unity  when  integrated 
over  the  snhere.  The  expectation  value  over  both  space  and  time 
of  the  amplitude  ~Kl  is 

„ p 

«^»  - z*  sl .c&pR  R1  bCR)  <>r.)  (38) 

Using  (31),  (35),  and  (37)  wo  obtain 

~ - n|  — j^x,  Co)  I5(a,  Si  +-  xt  (e)  X,  (r,s)j  (39) 

and 


A <a> 


1 1=>  / 


C iii  c=>> 


1\( 


o ®D 


^ . dJR*CcoS<S5) 

d <££}  S<A 


«d  «a>* 


i,  <».  *>  - >•  w 

I^(n,s)  = J J(k»R)  h(R)  )*  (>=.».  R) 

To  evaluate  J,  (n)  we  integrate  by  parts  and  obtain 

X,  I'O)  — - /a  CEO  cos  <E5D  R (cos  <S>) 

This  may  be  evaluated  (MacRobert.  1.948,  pg.  106)  to  give 


X,  (2  n)  = O j 


I ( r> . . \ 

i \ • j 


(40) 


rc^i)  (4i) 

r(pf o r ( °+ -z) 

To  evaluate  I,  Cr>;  we  note  that 

sLfi  (<os.<s&J  = ..  si„  ^ dfL£iss<a>) 

d 1:35  <=K  C os  '&i) 

which  vanishes  at  ©=oor  7r  , Thus  upon  integration  by  narts 


we  obtain 
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T 'rO  - - J (»>.)  (42) 

To  evaluate  X1(«J3!>  and  X.(o,s)  we  must,  assume  a form  for 
Vi  (R)  • Let  us  assume  equal  aprlori  probability  of  production  of 
loops  over  a sphere  of  radius  R.  ^ P and  zero  probability  outside. 
With  suitable  normalization  we  then  have 


KCR-)  = 


(■Zbz*  ^ <.  R. 

\ , 'H  7C  I \ „ 


1 


f A <z\ 

K O ) 


O for-  /R  > R. 

By  (41)  we  need  only  consider  mo de3  of  odd  order,  and  for  these, 

FR, 


~*I,C2-ljS)  jjM,  C uj  (44) 

To  evaluate  ia(2*+i.s)  we  integrate  by  parts  and  obtain 

I,(2^LS)  - 2 1,  (2**1,  S)  trj(2MljS)  (45) 

where 


‘‘•i*-.,  o. 


Consider  finally  (^r,4  l,s)  „ We  use  the  identity 

j-w)  - - -~-gL 

and  integrate  by  parts.  Repeating  the  process  n times  gives 


c 

j d-  j,..,  (v) 


\ 5 **  r> ) 2 ' u^_ 

4—  "(n-ofT! 


Noting  that  for  small 

>(^)  ' ' , (2^1)'!  = (*"+  >lz«-0  • • • 

we  have 

^2l5l*  T (?b.,^-2V  V"  2^!  (w  R.) 

3 1‘l2n‘1'-'  or 

Using  (42)  and  (45)  we  obtain 


(47) 


= - tri  (dcLi-^X.  (g?+  *)  ±jv(^l kill 

'(fc. 


(48) 
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Consider  now  the  expansion  of  B in  terms  of  ortnoronai 
poloidal  modes  _§ . Thdse  are  related  to  the  7T7,  by 

t7*  X,  = V »(v.  S ) = - y * S,  = ■>.'  S, 

From  (28)  we  compute  then 

c _ „ , ,i=  i.  Si-iai  <&">  (49) 

From  (48),  which  gives  > the  mean  value  of  c«,  , 

and  from  (4f 'we  find  that  the  expectation  value  of  B>-  is 

b-  = ^ r s.,.„ 

= - ^£c^3)J,(Zn-K,;fL  (4*  ftVsfc  fry  t ^ 1 . sj]  U..(k„...,R) 

" (k.<—>  fO*  Of.c.-.  r)  r; 

X (cc^  &>) 


where 

£ _ B.^VW^vtfa-P  (50) 

The  expectation  value  of  the  coefficient  FT  of  Ft!  C-‘os  <ss)  is 


, . , . 's  fax  L *)  i-Xs  , _ . » 

r^°-r— - £(4"*5)I'(*~';  & (51) 

To  compare  these  results  with  observation  we  compute  17 

and  17  , the  coefficients  of  the  dipole  and  octunole  moments  for  var- 
ious value  s of  R*.  . The  results  are  given  in  table  1 for  the 
magnitude  of  these  coefficients  at,  the  surface  of  the  core,  whence 
they  can  readily  be  converted  to  the  corresponding  values  at  the 
surface  of  the  earth. 


“25  - 


TABLE  i 

THE  COEFFICIENTS  OF  F?(cos<2>)  AND  F>(eo=,<S})  AT  THE  SURFACE  OF  THE 

CORE 


P 

rep.) 

e 

c 

0.5 

-0.0222 

-0.00086 

0.8 

-G.08IS 

-0,00235 

0.9 

-0.120 

-0.01125 

1.0 

-0.162 

+0.C251 

Nowvfc  have  shown  orevi ■uisl,T  that  the  loops  oroducea 
b,v  downgcing  currents  are  onDOsite  in  sign  to  the  loots  nrod'iced 
by  uogoing  currents*  we  rennrked  that  the  loops  nroduced  by  the 
downgoing  currents  were  formed  farther  from  the  surface  of  the  core 
than  the  loons  formed  by  rising  currents  and  therefore  would  make 
a much  smaller  contribution.  Let  us  now  show  quantitatively  that 
this  Is  the  case.  Let  us  assume  that  the  downgoing  loops  are 
wodi.ee  d only  out  to  Ff.  = P“/2  and  that  they  are  onposite  in 
sign  and  equal  in  number  to  the  loops  resulting  from  rising 
currents,  which  are  produced  in  the  region  Fy2<.R<P  . Table 

1 shows  that  7-.  ( P. .)  >>  TL,  (P/S)  for  R.  ==  R 

The  superposition  of  the  upgoing  and  downgoing  loops  gives 
a net  field  cf  17..,  (FP.)  — 7..  ('FV<2. 1.  The  ratio  of  the  dipole  to 
the  octupole  terms  at  R = P is  given  in  table  2 for  FP.  = o- 0!5 
ana  1-0  . We  expect  that  R1.  is  somewhat  less  than  P : It  is 
chysically  irrroossible  to  oroduce  a Iood  of  finite  size  exactly 
at  R = P*  ? the  toroidal  field  from  which  the  loops  are  formed 
goes  to  zero  at  R=P  (Bullard,  1954).  Observation  indicates 
a value  of  about  15.7  for  the  ratio  of  the  P component  to  the  R 
component  at  the  surface  of  the  core, 


(Elsasser,  1941)  which 
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agrees  with  our  calculations  for  F?-  Q . g, 57  p> 


TABLE  2. 

RATIO  OP  ZONAL  DIPOLE  TO  OCTUPOLE  MODES 
AT  THE  SURFACE  OF  THE  CORE 


There  is  probably  no  point  in  computing  higher  modes  than 
the  octupole  term  because  the  external  fields  of  these  mode3  are 
strongly  influenced  by  the  random  motions  near  the  surface  of  the 
corej  thus,  as  with  the  observed  quadruoole  field,  it  is  most 
likely  that  they  arise  mainly  from  random  fluctuations  rather  than 
from  the  average . feed  back  eft cct  discussed  here. 

In  the  above  numerical  calculations  we  have  neglected  the 
effects  of  a possible  inner  3olid  core  because  the  volume  of  such 
a sphere  is  small  compared  to  the  volume  of  the  entire  core,  and 
because,  as  has  been  shown,  Q Rf. ) decreases  so  rapidly  with 

FT.  . If  one  desired  to  take  into  account  an  inner  core,  it 
would  only  be  necessary  to  superpose  a third  set  of  loops  confined 
within  a radius  IP,  P/3  to  neutralize  all  other  loops  in  R<  F^.  . 
The  net  field  would  then  be  r,'mi  ( Ra)  — fl*..  (P*/  2)  f (■)  I7.«  (P, ) where  ^ 

Is  equal  to  . The  effect  of  ^ IZ..  (Ft,)  will 

be  small. 

Finally,  we  note  that,  with  a toroidal  field  which  vanishes 
at  the  equator  and  the  poles,  our  assumption  that  the  loops  are 
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6 . Reversal  of  Field 

Runcorn  (1954)  has  cited  geological  evidence  to  the  effect 
that  the  dipole  field  of  the  Earth  ha3  reversed  itself  in  the  past. 
The  question  naturally  arises  as  to  whether  such  a reversal  can  be 
incorporated  into  the  model  presented  in  this  paper.  Apparently 
there  is  a means  of  reversing  the  dipole  if  one  assumes  a sudder 
large  increase  in  the  convective  motions  within  the  core,  as  the 
following  qualitative  discussion  will  indicate. 

The  decay  of  the  magnetic  field  in  a dvnamo  mav  be  t hou  ght 
of  as  a slipping  of  the  lines  of  force  beck  through  the  fluid, 
equivalent  to  an  attempt  to  unwind  the  contortions  introduced  by 
the  motion  of  the  fluid.  In  a steadv  state  dynamo  it  is  physically 
obvious  that  the  slip  velocity  and  the  fluid  velocity  are  comparable, 
which  is  another  way  of  saying  that  the  characteristic  time  of  the 
circulation  of  the  fluid  and  the  decay  time  of  the  magnetic  field 
are  compa-able.  Indeed,  it  can  be  shown  that  the  regenerative 
operation  of  the  dynamo  requires  a certain  amount  of  phase  shift 
resulting  from  the  decay  ( Bcndi  and  Gold,  1950),.  the  decay  acts  in  a 
very  crude  sense  as  the  commutator  of  the  dynamo.  One  might  expect, 
then,  that  the  regenerative  process  of  the  dynamo  would  be  thrown 
out  of  gear  if,  for  instance,  the  characteristic  time  of  the  cir- 
culation of  the  fluid  were  to  be  suddenly  decreased. 

To  sec  just  what  will  go  wrong  in  our  dynamo,  1 et  us 
consider  a loop  of  flux  produced  near  the  surface  of  the  core  by  a 
rising  current  as  shown  in  figure  6.  For  normal  steady  operation 
of  the  dynamo,  we  find  that  during  the  time  that  the  loop  is  being 
produced  and  then  drawn  out  to  reinforce  the  toroidal  field,  it  will 
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diffuse  as  is  inaicated  by  the  broken  line  in  figure  6;  the  side  of 
the  loop  nearer  the  surface  of  the  core  will  diffuse  un  through  the 
surface  of  the  core  into  the  mantle  and  so  escape  the  nonuniform 
rotation  of  the  core.  Hence  it  will  not  contribute  to  the  toroidal 
field,  which  then  results  primarily  from  t he  side  of  the  loop  nearer 
the  center  of  the  core.  But  suppose,  or.  the  other  hand,  that  Just 
as  the  loop  is  being  formed,  the  velocity  of  the  fluid  increases 
in  order  of  magnitude.  The  diffusion  of  the  loop  will  not  have 
time  to  occur  and  both  sides  of  tho  loon,  not  Just  the  inner  side, 
will  be  subject  to  the  nonuniform  rotation  of  ths  core  and  contri- 
bute to  the  toroidal  field.  The  outer  side  of  the  'loon  is  in  a 
direction  such  that  it  will  degenerate  the  existing  toroidal  field 
near  the  surface  of  the  core.  This  degeneration  will  continue 
until  the  toroidal  field  reverses.  The  reversal  will  start  at 
the  surface  of  the  core  and  occupy  an  increasinglv  thick  shell j 
loops  produced  in  the  reversed  layers  will  degenerate  the  dipole 
field.  The  characteristic  decay  time  of  the  toroidal  field  is  of 
the  order  of  10 ^ years  for  6 =3*ios  mho/m  (Elsasser  1950)  j the 
characteristic  time  of  formation  of  a loop  is  600  years  for  a char- 
acteristic length  of  1000  km  and  a velocity  of  0.05  mm/sec.  Thus, 
we  see  that  the  reversed  toroidal  field,  if  once  begun  by  a brief 
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raid  persist  long  enough  to  completely 


degenerate  fend  ultimately  reverse  the  dipole  field 
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Appendix 

We  shall  treat  tho  case  of  a convective  stream  closed  in 
itself  but  so  that  the  motion  in  one  (say  the  upward)  direction  is 
strongly  concentrated  whereas  the  return  flow  is  spread  over  a wide 
region.  As  announced  in  the  text  shall  show  that  the  return  flow 
makes  oniy  a negligible  contribution  to  the  formation  of  magnetic 
looo  in  the  meridional  planes. 

We  use  the  local  cartesian  system  (^,  1,1)  and  assume  a 
toroidal  field  of  the  form  (15),  Let  the  velocity  field  of  the  fluid 
'£-£+-  *£  : u represents  the  convection  along  the  Jf  -axis 

together  with  the  associated  infiux  and  efflux  and  the  necessary 
return  flowj  wr  represents  the  rotation  of  the  fluid  about  the 
y -axis. 

Without  too  much  loss  of  generality  we  nay  take  = O 

and  set 

u,  = - u.ct)X(£>  v<i)  z'(«r;  , ^ = w.  (t;  x(t;  y'k;  zc$) 

where  the  crimes  denote  derivatives.  We  represent  \\#  by  w^.  = o and 

= --net)  i s(*)  j ^ = n (t;  # roo  sir) 

where  u.  (U  and  0.(0  are  representative  of  the  corresponding 
magnitudes,  provided  all  the  other  functions  are  suitably  normalised. 
The  forms  (3)  and  (i)  guarantee  that  V iS  and  V vanish. 

From  (l)  we  have  for  the  first-order  perturbation  of  the 
magnetic  field, 

/geo*  / 

The  second-order  perturbation  field  is 
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£(l)  = fj t 7*  [(*»  + 

This  field  consists  of  the  sum  of  four  terms  which  arise  aenending 
on  whether  the  first  stem  of  the  interaction  is  taken  to  involve 
either  l»  or  v*  , and  whether  the  second  ster  involves  either  j»  nr>« 
Two  of  these  terms  are  quadratic,  in  j alone  and  in  alone;  it 
is  physically  obvious  that  they  will  not  contribute  to  loops  in  the 
meridional  or  5^-olanes.  This  may  be  demonstrated  quantitatively 
by  noting  that  the  £ -comnonent  of  the  term  of  second  order  in  (.t) 
and  the  -component  of  the  term  of  second  order  in  fl  (t)  are  zero. 
Thus  the  terms  cannot  contribute  individually.  If  we  combine  the 
nonzero  J> -component  of  the  former  and  the  nonzero  | -comnonent  of 
the  latter,  we  obtain  a field  which  can  be  shown  to  have  no  net 
circulation  about  the  i^-axis.  We  therefore  omit  these  quadratic 
terms  and  keep  only  the  two  mixed  terms  which  involve  jj  In  one  step 
of  the  Interaction  and  in  the  other.  After  somewhat  lengthy 
nut  straightforward  calculations  this  mart  of  the  second-order  field, 
say  S (t  ),  Is  found  to  have  the  followinr  comnonents  in  the  meridional 
planes : 

fe,(U  - Ktw*)  (-  R(,W]  2\{)  5CP  BCp 

+ ([I^H  scs;  sew] 

' ( 

+ S(j)  Z(«BC«]j 

6,(0  - Ht) 

i lh> 

+ ^X(S)  (\J  - S X (5,  2(yj  SCSJ  BCs) 
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wh  ere 


f(t)  = fc )t/u„CO  “)  = / 

- • A • 


The  field  linos  of  and  S.  are  given  by 


A\  - ^ 


u. 


ej  § _ ci 

V,  W\ 


_ Jjr 


For  u»  we  obtain  the  family  of  curves 


(c) 


Y(0  ?CS)  = C, 


(d) 


For 


vw 


+ vl  - C, 


(e ) 


where  C,  and  C»  are  the  parameters  of  each  family.  The  differential 
equation  for  che  field  lines  of  is  not  readilv  integrated  except 
in  soecial  cases.  We  therefore  consider  the  case  of  a gaussian 
distribution 

x (g)=  ex^40  , h-) , z(.y)  - ^ J 

F?p;  - ( sent-  , B(jr;  = 

Then  (d)  becomes 

r\cxpf_  a>('V  + JO]  = c<  or  =a‘)hA  (T) 

L C 

and  (a)  and  (b)  become 


£*<o  - f(t)az[-jrica)f^(o7  -*p[-§  (fV^+yk;. 

Cl  J 


8,(t)  = KOB.  2 6 <oCV)  exp  [--*(£*+ r)  - £lZ±)‘l 


(g) 


wnere 


xcv,  -^(1  + e^)  !h) 

The  neutral  line  where  CO  and  £>s  ({,)  vanish  simultaneously 
will  be  denoted  by  i„(*v)  , Jf„  C >\)  • We  366  from  (g)  that 
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£ ,r  - - ~ r (n  ) - **  UJ  ( i) 

s.K,  - - » >-(\)  ~ T7X) 

Substituting  (g)  into  (c)  we  obtain  for  the  field  lines  of  the 

family  of  curves 


[y  - xc*o]2  + 


6~>  c y ) 

vr  ( v) 


a7  r 
S = Cj 


(k) 


These  are  ellipses  for  a1  with  centers  on  the  neutral  line.. 

For  ^ > a*  the  lines  form  hyperbolas. 

The  oroblem  of  defining  the  "strength**  of  the  magnetic 
loop  generated  offers  some  difficulty.  The  q -component  of  curl 
vanishes  on  integration  over  the  volume,  but  this  does  not  mean  that 
there  is  no  net  circulation  of  magnetic  flux  about  the  neutral  line. 
Perhans  the  simplest  way  to  demonstrate  that  such  a circulation  has 
been  generated  is  the  following:  We  compute  the  flux  across  the  p ^ - 

plane  above  X (\)  , and  show  that  at  least  some  of  it  bends  around  so 
as  to  penetrate  the  surface  X = X.  ( t ) , £ > O,  The  fluxes  are 


j dy  s. 

X.(\) 

Using  (g)  we  find 


§ -a 


= 


r 

'ij'di  &, 


hr; 


♦ OO  00 

cX,  = f(  {)£>'j\  [^(O.  exp  [-•*  C*C+  Xl)-?(S~b)‘l 

Ji(i>  J 

«>  UKV  exp  j-§  Wt  y;w;  . J (Sco. 

x ■§»  §‘) 


In  the  integrand  of  the  expression  for  , we  have  X ^ X.  (\) 

From  (i)  it  follows  that  the  integrand  is  always  positive.  In  the 
expression  for  the  sign  of  the  integrand  is  the  same  as  the 

sign  of  60 (r^)  . From  (h)  we  see  that  60(r^)  is  positive  if  r^1  < 
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otherwise  negative.  But  for  r^‘>al  the  gaussian  factor  is  very 
small,  sc  small  in  fact  that  the  integration  over  a‘  constitutes 

only  a few  percent  of  the  total  value.  Thus  the  integrals  are 
positive  and  and  have  the  same  signj  hence  there  is  a net 

circulation  of  flux  about  the  neutral  line  independent  of  the 
relative  magnitudes  of  a,  b,  and  c. 

Since  the  existence  of  a net  circulation  is  independent 
of  a,  b,  and  c we  shall  limit  the  evaluation  of  the  fluxes  to 
the  special  case  b = 0,  a = c.  Then  ^(n)  •=  O and 


cX>  = 


Jot" 

2Vs 


f CO  B.  a 
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